Preprint typeset in JHEP style. - HYPER VERSION 



NLO QCD corrections to tri-boson production 



T. Binoth a , G. Ossola\ C. G. Papadopoulos\ and R. Pittau c 

a The University of Edinburgh, School of Physics, Edinburgh EH9 3JZ, UK 

b Institute of Nuclear Physics, NCSR Demokritos, 15310 Athens, Greece 

Departamento de Fisica Teorica y del Cosmos, CAPFE, Universidad de Granada, E- 18071 
Granada, Spain 



Abstract: We present a calculation of the NLO QCD corrections for the production of three 
vector bosons at the LHC, namely ZZZ, W + W~Z, W+ZZ, and W + W~W+ production. The 
virtual corrections are computed using the recently proposed method of reduction at the integrand 
level (OPP reduction). Concerning the contributions coming from real emission we used the dipole 
subtraction to treat the soft and collinear divergences. We find that the QCD corrections for these 
electroweak processes are in the range between 70 and 100 percent. As such they have to be 
considered in experimental studies of triple vector boson production at the LHC. 
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1. Introduction 

For TeV collider physics hard multi-particle final states are ubiquitous and theoretical calcula- 
tions can not provide reliable predictions without taking into account higher order information. 
Unfortunately the evaluation of one-loop amplitudes with many external particles is technically 
very challenging, which motivated a priority list for one-loop computations relevant for the Large 
Hadron Collider at CERN, the so called Les Houches wish list Due to the relevance for LHC 
phenomenology many new avenues have been explored in the last few years, ranging from evaluation 
techniques of Feynman diagram || f|, ||, |J to unitarity based approaches in different variations 

Higher order QCD results have been provided recently for multi-boson production pp — ► 



ZZZ, WWZ, HHH processes M, 0, 0, M , in the context of weak boson fusion EO. 0. 0, 0, E4 



pp — > Hjj with effective gluon-Higgs couplings [25], gg — » Hqq [26], and pp — > ttj [27 
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In two recent papers |28|, 29], a new technique (OPP) has been introduced for the reduction 
of arbitrary one-loop sub-amplitudes at the integrand level [[KJ by exploiting numerically the set of 
kinematical equations for the integration momentum, that extend the quadruple, triple and double 
cuts used in the unitarity-cut method [[H], [32|, [33[]. The method requires a minimal information 
about the form of the one-loop (sub-)amplitude and therefore it is well suited for a numerical 
implementation. 

In the present work, the OPP reduction is applied to the calculation of the next-to-leading 
order QCD correction for the production of three vector bosons at the LHC This includes the case 
of ZZZ production, as well as the W + W~Z, W + ZZ, and W + W~W + production. The physics 
motivation for a reliable prediction of these processes is two-fold: firstly one is sensitive to quartic 
vector boson couplings and secondly the leptonic decays are prominent Standard Model backgrounds 
for multi-lepton and missing energy signatures present in many new physics scenarios. 

As the triple vector boson production is genuinely an electroweak process one can not expect 
that the inclusion of QCD effects leads to the reduced scale dependence typically seen in this 
kind of calculations. In contrary it can be qualitatively understood that the LO predictions show 
a relatively small sensitivity when varying the factorisation scale. This is because the parton 
distribution functions are called for x-values which are around the scaling region where one has a 
very mild Q 2 -dependence. After adding the order a s corrections one expects to observe a LO type 
scale variation in the added contribution. 

The production of three Z bosons has already been discussed by Lazopoulos et al. in Ref. [fTE] . 
We also presented some preliminary results in Refs.|], |34|. The W + W~Z case has been studied in 
Ref. ]17| for all combinations of leptonic final states. Results for W + ZZ and W + W~W + production 
have not been presented in the literature yet. 

Our calculation is composed of two main parts: the evaluation of virtual corrections, namely 
one-loop contributions obtained adding a virtual particle to the tree-level diagrams, and corrections 
from the real emission of one additional massless particles from initial and final states, needed 
in order to control and cancel infrared singularities. The virtual corrections are computed using 
the OPP reduction [28, In particular, we make use of CutTools [3~5|, a FORTRAN90 code 
that implements the general method of reduction. Concerning the contributions coming from real 



emission we used the dipole subtraction method |36| to isolate the soft and collinear divergences 
and checked the results using the phase space slicing method p7L |38| with soft and collinear cutoffs, 



as outlined in [39, 4C 



The paper is organized as follows. In Section |2|, we report the details of the calculation of the 
virtual part. Section |3| is devoted to the discussion of soft and collinear singularities. In Section [|, 
we show our results, including transverse momentum and rapidity distributions for the different 
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processes studied in this paper. Finally, in Section [^, we will give a summary of the work done and 
present our conclusions. 



2. Virtual corrections 



We consider the process 

q(pi) + qfo) — ► v( P3 ) + v(p A ) + v( P5 ) 

where V = Z, W. All momenta are chosen to be incoming, such that YliPi = 0- 



(2.1) 




Figure 1: Tree-level structures of Feynman diagrams contributing to qq — > VVV, where V = Z,W. 
Dashed internal lines can represent W, Z, Goldstone bosons or photons. 

At the tree-level, diagrams can be grouped in four different topologies, which are illustrated in 
Fig. [I]. One-loop corrections are obtained by adding a virtual gluon to the tree-level structures, as 
depicted in Figs. |2] and [| 

We perform a reduction to scalar integrals using the OPP reduction method [^, ^9j. In this 
approach, we need to provide the numerical value of the numerator of the integrand in the loop 
integrals. We refer to it as the numerator function N(q), where q is the integration momentum. 

The numerator function N(q) can be expressed in terms of 4-dimensional denominators Di = 
(q + pi) 2 — mf as follows 



m—1 



N{q) = ^2 [d(«o«i^3) + d(q; z zi 

«0<*1<*2<*3 
m—1 

+ [cO*oV*2) +c(q;i i 1 i 2 )] \ [ D 

«0<«1<«2 i^t 
m—1 m—1 

+ [ b (ioh) +6(g;«o*i)] II D 



m—1 

'2':s)| II Di 

m—1 



io<h 
m—1 



m—1 



+ ^2 [ a (^o) + «(<?; «o)] Yl Di ■ 



(2.2) 
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The quantities diiQi^iz) are the coefficients of 4-point scalar functions with denominators labeled 
by i , ii, i 2 , and i 3 . In the same way, c(i ^i^2), K^oh), and a(i ) are the coefficients of the 3-point, 
2-point and 1-point scalar functions, respectively. The other quantities appearing in Eq. ( |2.2| ), 
marked with a "tilde", vanish upon integration over q. Such a separation is always possible and 
the set of coefficients d, c, b, a is immediately interpretable as the ensemble of the coefficients of all 
possible 4, 3, 2, 1-point one-loop functions contributing to the amplitude. 

Since the structure of Eq. Q2.2p is general, namely independent from the particular process that 
we want to study, the task of computing the one-loop amplitude is then reduced to the algebraical 
problem of fitting the coefficients d, c, b, a by evaluating the function N(q) a sufficient number of 
times, at different values of q, and then inverting the system. That can be achieved quite efficiently 
by singling out particular choices of q such that, systematically, 4, 3, 2 or 1 among all possible 
denominators Di vanishes. Then the system of equations is solved iteratively 1 . 

First one determines all possible 4-point functions, then the 3-point functions and so on. In 
summary, simply by evaluating the numerator function N(q) for a given set of values of q, we can 
extract all the coefficients in Eq. ( |2.2|) . 

As a possible future development, the numerical evaluation of N(q) could be performed auto- 
matically without relying on Feynman diagrams, by means of recursion relations. For the current 
project however, we still follow the traditional approach of computing all the expression originat- 
ing from Feynman diagrams, and use them to evaluate numerically the numerator functions. An 
example in Section |2J] will clarify the details of the technique used. 

The coefficients determined in this manner should be multiplied by the corresponding scalar 
integrals. Since, in the process that we are studying, no g-dependent massive propagators appear, 
we will only need massless scalar integrals. They are computed using the package OneLOop written 
by A. van Hameren ||. 



The last step is the calculation of the rational terms. As explained in Ref. f42 |, there are two 
sources of the rational terms: a first contribution, that we call R±, originates from considering the 
fact that the denominators appearing in the scalar integrals are n-dimensional objects, while the 
expansion of Eq. ( |2.2|) is purely 4-dimensional. These contributions can be automatically extracted 
in the reduction process, either by computing extra-integrals as explained in Ref. [J29|| , or by means 
of a modified version of Eq. (|2.2| ) in which the numerator function is expressed directly in terms of 
n-dimensional denominators. The second approach is illustrated in Ref. ^2| and implemented in 
the package CutTools [35|. We checked that the results obtained for Ri with the two methods are 
in perfect agreement. 

The second contribution, that we call R2, is instead originating from the numerator function. 



*A method to optimize the solution of the system has been very recently presented in Elf 
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For many processes, N(q) can be treated as a purely four-dimensional object. However, in general, 
it should be written as N(q) = N(q) + N(q 2 , e), where N(q) is the n-dimensional numerator function. 
N(q 2 , e) can originate, for example, from the contraction of Dirac matrices or from powers of q 2 in 
the numerator function [(43| and vanishes in the e — ^ limit. In Ref. [42 we discussed in detail this 
topic and showed how R2 can be obtained, for example, by using a set of tree-level like Feynman 
rules. For the calculation presented in this paper, however, it is easy to extract these remaining 
rational parts directly. 

2.1 ZZZ production 

In this subsection we describe the evaluation of the virtual QCD corrections to the process qq — » 
ZZZ. The virtual part of the calculation involves eight different diagrams, which have been depicted 
in Fig. 0. Each diagram should be evaluated for six permutations of the final particles. 





Figure 2: Diagrams contributing to virtual QCD corrections to 



ZZZ 



As an example, let us consider the pentagon diagram (the last diagram of Fig. [|). In our 
notation, the integrand will read 

N 5 (q) 



with 



Ab 9 [q 2 ] [(<? + Pi) 2 } [{q + Pi+ P5) 2 } [(q -P2- P3) 2 } [(q - P2) 2 } 

N 5 (q) = {u( P2 ) 7 Q P {q _ p2) V z z P {q . P2 . pz) V? P (q+Pl+P5) Vi P (q+pi) r uipx)} 



(2.3) 
(2.4) 



The function P(q) is the numerator of the quark propagator 

P (g) = i > 

while = V z ■ €i , namely the contraction between the polarization vector of the i-th Z boson £j 
and the 7-matrix in the vertex Zqq 

V* = ie^{gju- + gju+) (2.5) 
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where 

s _ s 2 Qf 
gf = — Q f , gj = , cj ± = (l±7 5 )/2, s = sin^v, c = cos6 w . (2.6) 

For any fixed value go of integration momentum, and for a given phase space point, N 5 (q ) is simply 
the trace of a string of known matrices. After choosing a representation for Dirac matrices and 
spinors, we evaluate N(q) by performing a naive matrix multiplication. By providing this input to 
the reduction algorithm, we can compute all the coefficients of the scalar integrals (in other words, 
the "cut-constructible" part of the calculation). 

In the same fashion, we can repeat the calculation for the other seven diagrams. However, our 
method allows for a further simplification: for each fixed permutation of the final legs, only the 
q-dependent denominators of Eq. (|2.3|) will appear also in the remaining diagrams. Therefore, we 
can combine all diagrams in a single numerator function and perform the reduction directly for the 
sum of such diagrams, allowing for a one-shot evaluation of the resulting scalar coefficients. 

We checked that our results, both for poles and finite parts, agree with the results obtained by 
the authors of Ref . . 

2.2 W+W~Z production 

With the same technique we also evaluated the virtual QCD corrections to the process qq — > 
W + W~ Z. The structure of the diagrams is more involved with respect to the ZZZ case. There are 
in fact 19 different tree level diagrams. Adding QCD corrections, we obtain 58 one-loop diagrams 
contributing to this process. In addition to the structures already depicted in Fig. ||, in this case 
we can also have the diagrams appearing in Fig. |3|. 

A very similar calculation has been presented recently by Hankele and Zeppenfeld ]17|]. They 
studied the NLO QCD corrections to the production of 6 leptons in hadronic collisions, via WW Z 
production. A comparison with their results, however, is not straightforward and has not been 
performed yet. 

2.3 W + ZZ production 

Concerning the production of W + ZZ, we have 15 tree-level diagrams, which, after adding QCD 
corrections, give rise to 69 diagrams at NLO. Since in this process we generate a single charged 
W + , the initial state should be of the type ud (rather then uu as for W + W~ Z and ZZZ). 

2.4 W + W~W + production 

Starting again from an initial state of the type ud, we should consider 15 diagrams at the tree-level 
and 53 diagrams including NLO QCD corrections. 
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Figure 3: Additional NLO structures contributing in the production of W + W~ Z , W + ZZ, and 
W + W~W + , that do not appear in the ZZZ case. Dashed internal lines can represent W, Z, Goldstone 
bosons or photons. 

3. Real emission 

The real emission corrections for the production process of three vector bosons 

q + q^V + V + V (3.1) 

fall in the following three categories 

q + q^V + V + V + g (3.2) 

g + q^V + V + V + q (3.3) 

g + q + V + V + q . (3.4) 

IR divergences arise if a massless final state particle becomes soft or collinear to an initial parton. 
We deal with the IR part of the calculation by using the two cut-off phase space slicing method 
|39| , [40ll and the dipole formalism of Catani and Seymour [|36|. Let us first provide the relevant 



formulas for the dipole subtraction method. 



3.1 Dipole subtraction 

The partonic cross section at the NLO level consists of Born term (B), virtual corrections (V), 
collinear counter terms (C) defined on the 3-particle phase space and the real emission corrections 
(R). Dipole terms (A) which approximate the real emission matrix elements in all soft/collinear 
regions are subtracted from the real matrix element before integration over the four-particle phase 
space. The same terms are added back, integrated over the dimensionally regulated phase space of 
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the soft /collinear particle: 

I [da*-daf^ + J [da*+daV + Jdrf 9 + da°] . (3.5) 



„NLO 



VVVg VVV 



After subtracting the dipole terms the real emission cross section is finite and can be evaluated in 
4 dimensions. The same is true for the other terms after the pole parts have been canceled. 
The colour averaged leading order contribution is given by 



d <=jf^-\ MB \ 2d ®wv ( 3 - 6 ) 



where A4 B is the kinematic part of the leading order amplitude and s i2 = 2pi ■ p 2 . If two (three) 
vector bosons are identical a symmetry factor C s = 1/2 (C s = 1/6) has to be included. The three 
particle phase space of the vector bosons is denoted as d& V vv- The real emission corrections are 
defined on the four particle phase space d^vvvx where x can be either g, q, or q. 

In the case of a qq initial state two dipoles are needed as subtraction terms. The subtraction 
term for the gluon emission off the quark (neglecting 0(e) terms) is 

where 

pi • p 2 - p 2 ■ Pe - Pi ■ Pe 

x 

Pi ■ P2 

Pie = xpi , K = p 1 +p 2 -p 6 , K = p 16 +p 2 

9 (K + ky 

Pj = A Pj (3.8) 

defines the dipole kinematics: q(pw) + q{p2) — > V{p?) + V{p^) + V{p$). The subtraction term for 
gluon emission off the anti-quark is obtained by interchanging p 1 <-> p 2 . The real emission cross 
section including subtraction terms reads 



do* -da* ~ ( 



qq qq N 2s 



12 



C F \Mfg\ 2 - V qi96 > q2 - £>^6,<?ij d^ VVVg (3.9) 



The part of the NLO cross section which is defined on the 2-^3 phase space is obtained after 
analytic integration of the dipole terms over the phase space of the unresolved particle. A collinear 
counter term is added to treat the collinear 1/e pole which is absorbed into the parton distribution 
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functions at a scale jj F . All details can be found in [pq| . The part which has to be added to the 
virtual corrections is given by 



2 3 2tt 



2-, 



' ' J ' f dx }C q > q (x) dafJx Pl ,p,) + ^ / dx m{x) do* fa, xp 2 ) 



2tt 



2tt 



(3.10) 



where the term 

/C M (x) = /C M (x) 



1 + x 2 
1 — x 



4 log(l — x) 
1 — x 



1 - x) - 2(1 + x) log(l - x) 



contains plus distributions which are defined as usual 



i 

/ 



dx 



1 — X 



f(x) = / dxg(x) 



,/( *)-/(!) 
1 — X 





For initial states with a gluon no soft contribution is present and thus one has 



a 



NLO 
91 



+ 



~ dof q 



VVV q VVVg 

In this case only one subtraction term is needed, namely 

C s 1 
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T R \M^,\ 2 - V 9iq6 ' g2 



gq 1 



VVVq , 



where the dipole is given by 



T> 9iq6 ' q2 = ?™' Tr [1-2x(1-x)] \M*(p 3 )\ 2 



x 2pi -p e 



(3.11) 



(3.12) 



(3.13) 



(3.14) 



The momentum mappings pj are identical to the ones in Eq. 

The initial state collinear singularity is again absorbed by the pdfs through a counter term 



i 

/ K = f ^w^do^xp^) 







Sl2 



IC 9 ' q (x) = [x 2 + (1 - x) 2 ] log ( ) + 2x(l-x) + 2[x 2 + (l-x) 2 ]log(l-x) . (3.15) 
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The formulas for the cases qg, qg, gq are identical up to relabeling of momenta. 

The hadronic differential cross section with hadron momenta Pi and P 2 is the sum over all 
partonic initial states convoluted with the parton distribution functions 

d(T (Pj, P 2 ) = ^ J dz 1 dz 2 f a (zi,HF)fb(z2,HF)d(7ab(ziPl,Z2P2) (3.16) 
ab 

The sum runs over the six partonic configurations qq, qq, gq, qg, gq, qg. 
3.2 Phase space slicing 

To have an independent check for the real radiation we have also implemented the phase space 



slicing method in its two cut-off variant |39], f40| . One splits the phase space in soft, collinear and 
hard regions with the help of the cut-off parameters 5 S and S c . In the soft region the 2 — > 4 matrix 
element is replaced by the eikonal approximation. In the collinear region one has a convolution of 
a splitting function with the Born term. Adding the soft /collinear parts to the virtual corrections 
all poles cancel and one obtains the three-particle contribution 



a® 



S / dZ ^ d(T ab [fa(zi,» F )Mz 2 , lip) (A S Q + A% + 2^ C ) 
n a,b J 



with 



+ f a (zi,^ F )f b (z2, IX F ) + fa(Zl,IX F )fb(z2, fX F ) (3.17) 



A s = A\n 2 6 s C F 

A s c (q - qg) = C F (2 \nS s + 3/2) In ^ 

fx f 

K = ■ (3.18) 

aa ab 



The / functions are given by 

" l - 5 ^ dz 



fa(x,fX F ) = ^2 —f b (x/z,fX F )Pab(z). (3.19) 

L J X Z 



Pab(z) = Pab(z) In ( 5 C X — 2xpi 2 P2 ) _ P > ab ( z ) . (3.20) 



where 



Z Ll F 

The upper limit 1 — 5 S ensures that the soft region which is already dealt with is excluded. The 
Kronecker 5 a b indicates that for a ^ b there is only a collinear divergence and no soft cut-off is 
needed. 
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In our case we need only the splitting functions P qq (z) and P gq (z). If we write P ab (z,e) = 
Pab(z) + eP' ab {z), we have 

P gq (z) = C F ^-^ (3.21) 
P' qq {z) = -C F (1 - z) 
P gq (z)=T R (z 2 + (l-z) 2 ) 
P' gq {z) = -2T R z(l-z) 

We see that the / functions contain an explicit logarithm of S c as well as logarithmic dependencies 
on 5 S which are built up by the integration on z\, z 2 . 
The four-body contribution is given by 



a,b=q,q,g 

with the hard-non-collinear partonic cross section given by 



^2 / dz 1 dz 2 f a (z 1 ,fi F )f b (z2,fi F )da^ b , (3.22) 

7Z <J 



da ab = [ _ Y^\M ah \ 2 d<$> VVVx} (3.23) 

^S\2 JHC 

where Xll-^^l 2 ^ s ^ e two-to-four body squared matrix element averaged (summed) over initial 
(final) degrees of freedom, d^vvvx is the four-body phase space and the hard non-collinear region 
denoted by HC is defined by 



E 6 > 5 S 1— 

2p 1 ■ p 6 , 2p 2 ■ p 6 > S c S12 (3.24) 

where p§ is the momentum of the soft /collinear parton with energy £V 

Both methods have been implemented in HELAC(45|]. The results show excellent agreement 
between the two methods. In the numerical results presented below we only show the results of the 
dipole subtraction approach. 

4. Numerical results 

We present in this Section a selection of the results that we obtained for the four processes studied 
in this paper. 

The complete virtual part of the next-to-leading order calculation for the four processes studied 



in this paper has been performed using CutTools [35] and also checked against an independent code 
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The two programs provide identical results for the amplitudes studied. As further tests, we checked 
that the tree-level results obtained using Feynman diagrams coincides with the results obtained with 



HELAC |^5| and that we reconstruct the correct structure for the poles after integration. Concerning 



the finite parts, we agree with the results obtained by the authors of Ref . |y| , for the production 



of three Z bosons. In this section we will mostly focus on the processes for which no results have 
appeared yet in the literature, namely the production of W + W~W + and W + ZZ. 
We use the following values for the electroweak parameters: 

M w = 80.4 GeV , M z = 91.1875 GeV , G F = 1.16639 • 1(T 5 GeV~ 2 . (4.1) 

In all cases presented here, we set y/s = 14 TeV and used CTEQ6L1 [44 ] with a s (Mz) = 0.129 at 



NLO. For the electroweak couplings we use the the G M scheme with 

v^M^sin 2 ^ 



71 

and 

2 



sin' 6 W = 1 - M^/M| (4.3) 



The tree- level cross section has been evaluated using the HELAC event generator ^5[. The same 
programme, appropriately adapted, has been used to calculate also the real corrections. The virtual 
corrections have been calculated on the basis of unweighted tree-order events produced by HELAC 
with an indicative CPU time of 180 ms per event, which is quite good taking into account that the 
numerical calculation of one-loop amplitudes (the numerators of the OPP method) is performed 
using standard momentum representation of Feynman graphs without any optimization. A conser- 
vative comparison with the efficiency of the tree order calculation, based on HELAC, shows that a 
further improvement of the order of 10 1 — 10 2 is to be expected. 

Since the purpose of our paper is to show the feasibility of the OPP method in a realistic 
situation, the results we present here are indicative and they constitute by no means a detailed 
discussion of the phenomenology of these processes. Partial results have been already presented in 
ill- 



It should be mentioned however that all results are available as (un)weighted events, which 
means that an exhaustive study in the full phase space, both for three- and four-particle final 
states 2 , poses no problem and will be postponed for the future, taking into account also decay 
products and intermediate Higgs contributions. 

In Figure (|, we show results for the pt distributions of all processes. For each phase space point, 
the pt of each of the three bosons gives an entry in the histograms. The final result is then divided 



2 Of course both positive and negative contributions have been taken into account, separately. 



12 



Process 


scale \i 


Born cross section [fb] 


NLO cross section [fb] 


zzz 


3M Z 


9.7(1) 


15.3(1) 


wzz 


2M Z + M w 


20.2(1) 


40.4(2) 


wwz 


M z + 2M W 


96.8(6) 


181.7(8) 


www 


3M W 


82.5(5) 


146.2(6) 



Table 1: Cross section for the four processes, corresponding to the distributions in Fig ||. Different values 
of the factorization (renormalization) scale are used for the different processes. 

by 3, yielding, as a normalization factor, the total cross section. The corresponding K- factors are 
depicted in Figure |5|. In the W + ZZ and W + W~W + cases, we observe an interesting increase in 
the K-factor for high values of the transverse momentum. 

The corresponding total cross sections are contained in Table [I]. As we can see the NLO 
corrections are quite significant, resulting to overall K-factors of order ~ 2. 

In Figure ||, we show, as an indicative case, the rapidity distribution for WWW production, 
which is the process with the highest cross section. Also here, each of the three bosons gives an 
entry in the histograms, that are eventually normalized to the total cross section. 

The K-factor appears to have now an important dependence on the phase space, especially near 
the borders of the available rapidity region 

Let us discuss now the results obtained for the production of W + W~W + and W + ZZ . In 
Tables |2] and |3] we present the results for the cross sections (in fb) of pp W + W~W + and pp — > 
W + ZZ, respectively. Each table contains the Born level, the NLO result and the corresponding 
K-factor. 



scale 




&NLO 


K 


H = M/2 


82.7(5) 


153.2(6) 


1.85 


/i = M 


81.4(5) 


144.5(6) 


1.77 


/i = 2M 


81.8(5) 


139.1(6) 


1.70 



Table 2: Cross section pp — > W + W W + in fb for different values of the factorization (renormalization) 
scale. In the table above we set M = 3M Z - 



5. Summary and Conclusions 

In this paper we considered the production of three vector bosons at the LHC We discussed four 
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Figure 4: Transverse momentum distribution, as denned in the text, for the four processes pp 
NLO (solid line) compared with the LO contribution (dashed line). 



VVV: 



processes, namely ZZZ, W + W zT, W + ZZ, and W + W W + production: for each process we cal- 
culated the next-to-leading order QCD corrections, presenting our results in the form of transverse 
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momentum and rapidity distributions. The QCD corrections are quite sizable, with a K-factor of 
order 2. The K-factor is rather uniform in distributions, while shows an important dependence 
on the phase space as far as rapidity distributions are concerned. Given the size of the corrections 
the QCD corrections have to be taken into account in experimental studies at the LHC. 

This paper also represents the first complete calculation of physical cross-sections performed 
using the recently introduced OPP method for the reduction of one-loop amplitudes, in which the 
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y y 

Figure 6: Rapidity distribution, as denned in the text, for pp — > W + W~W + : on the left plot, NLO (solid 
line) compared with the LO contribution (dashed line) in logarithmic scale; on the right, the corresponding 
K— factor. The scale is set to fi = 3M\y. 



scale 




&NLO 


K 


/i = M/2 


20.2(1) 


43.0(2) 


2.12 


fi = M 


20.0(1) 


39.7(2) 


1.99 


H = 2M 


19.7(1) 


37.8(2) 


1.91 



Table 3: Cross section pp — > W + ZZ in fb for different values of the factorization(renormalization) scale. 
In the table above we set M = 3Mz- 

reduction to known integrals is performed at the integrand level, using the Fortran code CutTools. 

The efficiency of the OPP method is quite good. It can be further improved by developing the 
numerical evaluation of the integrand in the one- loop amplitude by means of recursion relations [|46| . 
without relying on Feynman diagrams. 

We conclude that the OPP method is a viable alternative to perform phenomenologically rel- 
evant one-loop calculations, as it does not rely on the recursive evaluation of scalar and tensor 
momentum integrals. Its versatility and simplicity make it a very good candidate for the construc- 
tion of a universal NLO calculator/event-generator. 
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